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Abstract 

In this work we construct a general class of exactly solvable non- 
relativistic bi-dimensional quantum systems with position-dependent 
masses (PDM). These systems are isospectral to a given system with 
constant mass. The case of a charged particle with a PDM interacting 
with an external magnetic field is included in the present investiga- 
tion. We apply the approach in order to construct the SU(2) coherent 
states in some examples which are isospectral to the two-dimensional 
anisotropic harmonic oscillator, and discuss the impact of the intro- 
duction of special non-homogeneous external magnetic fields. 



*E-mail: dutra@feg.unesp.br 



1 



1 Introduction 



The interest in the problem of position-dependent mass (PDM) systems has 
been growing in the last years, both from the non-relativistic and the rel- 
ativistic point of view [1]-[3T]. In fact, there are many applications to dif- 
ferent problems of physics like quantum dots [32l |33], compositional graded 
crystals |3l], quantum liquids [35], metal clusters [36], neutron stars [37] . 
among others. In particular, recently the Wigner function for some classes of 
position-dependent Schroedinger equations were constructed and analyzed in 
the one-dimensional case [1]. The approach used in that work is capable to 
generate a class of position-dependent mass systems, which is isospectral to 
a given exactly solvable potential with constant mass. In this work, we will 
extend that approach to higher dimensions. Moreover, as far as we know, 
the majority part of the works dedicated to the research of PDM systems [1], 
deals with one- dimensional systems. Although, there are physical systems 
like those where a magnetic field [SS] is present, which leads naturally to the 
need of a two-dimensional analysis. Finally, to our knowledge, no work in 
this subject has discussed the case of PDM in the presence of magnetic fields. 
In this work, we intend to partially fill this gap. 

The method will be applied to the case of the anisotropic two-dimensional 
harmonic oscillator as well as in the case of those systems under the influence 
of some non-homogenous external magnetic fields. The SU{2) coherent states 
will be constructed for a number of systems which we choose to illustrate our 
results. 

In order to get the exact solutions for the two-dimensional Schroedinger 
equation we trace two routes. We begin by performing a general spatial 
variables change in the case of systems which are not under the effect of 
magnetic fields. Then, in the case of magnetic interaction, we begin by 
performing a time-dependent variable transformation followed by the spatial 
transformation. In the first case, we consider three examples of real variable 
transformations, which are respectively: the so-called polynomial one, the 
one using elliptic cylindrical coordinates and, finally, the bipolar coordinates 
transformation. Once we have the exact solutions for the eigenstates of 
the Schroedinger equation, we proceed with the construction of the SU{2) 
coherent states [H [391 SOI SI] • These last allow us to acquire a notion of the 
respective classical behavior of those position-dependent massive particles, 
similarly to what happens with the Wigner functions in one spatial dimension 
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This work is organized as follows: In the section 2 we present the approach 

wc arc going to use. Then, in the section 3 we apply it to some special mass 
dependencies, constructing classes of two-dimensional PDM systems which 
are isospectral to the anisotropic bi-dimensional harmonic oscillator. In the 
section 4, the case of magnetic field is taken into account. Finally we trace 
our final comments in the section 5. 



2 Class of isospectral two-dimensional position- 
dependent mass quantum systems 



In this section we present the approach which is capable to generate a class 
of models with position-dependent masses from a constant one, and which 
can also include the interaction with a magnetic field. In this case the two- 
dimensional Schroedinger equation is given by 



where the constant mo is the mass of the particle in the original system, x 
and y are the corresponding spatial coordinates, Px and Py are the respective 
momenta. Furthermore, we will work with magnetic fields in the Coulomb 
gauge, where — 0. Now, performing a general variable transformation 

in the spatial coordinates as 



2 mo 




(1) 




(2) 



one gets 



u 



fv dx + gv di 



(3) 



which can be written as 






Jacobian matrix, and its Jacobian is J = 



fv d , fu d 



After straightforward calculations, one can obtain the following expression 
for the Laplace operator in the transformed coordinates, 



xy 



91 +9'y = Y2 [{al + fl) dl + {gl + fl) dl - 2 {g^g, + + 



- f d 



gv du ) - gu dv (^y ) + du 

gu dy (^y ) - gv du (^y ) + fudv l^] - fvdu 



k 
J 

fu 



k 
J 

fu 

T 



du + 
du. (6) 



In order to achieve a system where we have a new Schroedinger-type equation 
with a position dependent mass one must require that the crossed derivative 
term 5^^, must vanish, leading to the restriction 



gugv + fu fv = 0. 
Afterwards, we should also impose that 

gl + fl = gl + fl 



(7) 



(8) 



to guarantee that the mass term is the same in the both terms and d^. 
Note that the equation ([71) once solved leads to 



fu 



gu gv 

' fv 



(9) 



provided that 7^ 0. By substituting dHD in ([H]), one gets 



2 I ^2 2 I 2 I 9v 
9u + fv = 9u + 9u{j 



9t (1 + 1 



gl 



fl + gl 
fl 



Thus, from (fTO!) we conclude that 

fv~gui fv — ^gu- 

Substituting (fTT!) in ([9]), one obtains 

gu gv 



fu 



±gu 



Tgv 



(10) 



(11) 



(12) 
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In fact, in all the cases considered here, the terms linear in the derivatives 
in the transformed Laplace operator (jSj) vanish. Moreover, the transforma- 
tion will keep = for all the cases considered. Thus the Schroedinerer 
equation in those transformed variables is written as 



M{u,v) 



ip {u,v) = E ip {u, v) 
(13) 



J2 



9l + fl 



Vif,g) + 



2 M{u,v 



-A\f,9), (14) 



2 M{u,v, 
where 

M {u, v) = rriQ 

with A = M (m, f ) ~t {f{u,v),g{u,v)). 

As an illustrative example, we will first deal with two-dimensional position- 
dependent systems which are isospectral to the anisotropic harmonic oscilla- 
tor and, then, with the ones which are isospectral to the case of a isotropic 
harmonic oscillator under the influence of homogeneous magnetic fields. 



3 2D PDM systems isospectral to the anisotropic 
oscillator 

As advertised, the first example which we will present here is the one of the 
anisotropic oscillator with constant mass governed by the equation 

- TT^ V V + {uf + ul y')^ = E ^. (15) 

In this case the eigenfunctions can be straightforwardly obtained, and are 
given by 



1 



^2(m+n+i) Trnlml X Y 



Hr, 



y2/(«,' 



X 



Hr, 



y/2g{u,v) 



Y 



X 



xexp 



f{u,v)Y fg{u,v) 



X 



Y 



(16) 



where X = ^/ihj (mo Wi), Y = ^/2h/ (mo ^2), UJi = q and UJ2 = P- 
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At this point we could make a study of the eigenfunctions and eigenvalues 
for the PDM systems which are isospectral to this one. However, we pre- 
fer to construct the so-called SU{2) coherent states, as defined in |39]-|^. 
Those states present the interesting feature of having their highest proba- 
bility density over a trajectory which corresponds to the classical one when 
7i — 7- 0. Furthermore, as the eigenstates, they are stationary wave-functions 
in contrast with the usual coherent states. As we are going to see, they 
will lead us to the conclusion that the isospectral PDM states will present a 
behavior which is very similar to the one of their "parent" anisotropic har- 
monic oscillator. Those coherent states can be written by using the definition 
introduced by Chen et al. [3S]-|1I] 



1 ^ f L 

<^{u,v,t) = - — -^V ^ r^ipnm{fiu,v),g{u,v)), (17) 

(1 + rr) ^ K=o^^ 

where the quantum numbers are defined as n = p K, m = q{L — K), with 
K = 0,1, 2, L, p and q being integer numbers. The complex parameter is 
such that T = Ae^'^, where = |, written in terms of polar coordinates, is 
used in order to make the connection with the classical trajectory. 

From now on, we will devote this section to develop explicit examples. 
As our first example of transformation functions, we deal with the second 
degree polynomials like 



11 11 

f{u,v) = -aiu'^+-biv'^+ciuv+di, g{u,v) = -a2U^+-b2v'^+C2UV+d2. (18) 

After imposing the restrictions ([7]) and (8), we are led to 



f{u,v) = T^C2u'^±^C2V^ + CiUV + di, g{u,v) = ±^Ciu'^T^CiV^ + C2uv + d2, 

(19) 

where one can see that only four arbitrary constants are left, and one have 
that 

fu+9l = f.+9l = {ci + cl)iu^ + v% (20) 
and the Laplace operator becomes 
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vL = -j^ M + fl) dl + {gl + fl) dl] 



(21) 



from which one conclude that the position-dependent mass is given by 

M{u,v)^mQ{cl + cl){u +v^). (22) 
In this case, the effective anisotropic potential looks like 

2 



V{u,v) 



mo 



(u^ — v'^) +CiU V + di^ uf + (v^ — U^) + C2U V + 



(23) 



and the wavefunctions are written as 



■■H„ 



V2{ 



£2 (,,2 



U — V 



ciu V + dij 



Pi 



X 



Hr, 



exp 



exp 



V2 (f {V^ - K^) + C2UV + d^) 



[f (w^ - v'^) + C1UV + di) 
Pi 

(f {V^ - U^) + C2UV+ ^2) 



X 



X 



(24) 
(25) 



In the Figures 1 and 2 (we use A = h = mo = 1 and L = 20, when generating 
the plots presented throughout this work) , one can see the behavior of the sys- 
tem in some typical situations. At this point it is interesting to note that the 
transformation functions can be chosen as the ones related to the parabolic 
cylindrical coordinates in a plane, where f = u v and g = {u^ — v'^) /2, which 
generate the following Laplace operator 



* uv 



1 



(u + 



(26) 



This make us remember that one could work with a general class of solutions 
as the ones coming from the orthogonal coordinate systems, where the metric 
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is diagonal. As examples we will consider the case of the elliptic cylindrical 
coordinates on a plane where 

f{u,v) = a sinh(-u) sin(f), g{u,v) = a cosh.{u) cos{v). (27) 
In this case, the Laplace operator is written as 

2 



(28) 



a?[cos{2v) — cosh(2-u)] 

with the mass M {u,v) = mo a^[cos(2v) — cosh(2M)]/2. The corresponding 
effective potential is now given by 



V{u,v) 



a mo 



[sinh (u^) sin (v) ul + cosh {u^) cos {v) ul\ 



(29) 



The wavefuntions are written in this case as 



V'nm {U, V) 



--H„ 



\/2[a sinh(M) sin(t')) 



y2("»+"-i)m!n!PiP2 

-\/2(a cosh(-u) cos(t')) 



Pi 



X 



exp 



P2 

a cosh(M) cos(v) 



cxp 



a sinh(M) sin(t') 
Pi 



38) 



whose typical behavior appears in the Figure 3. 

Now, in the bispherical coordinates on a plane, the transformations are 



a sinh(u) a sin(^;) 



cosh(ii) — cos(v) ' cosh(w) — cos(i') ' 

and the corresponding Laplace operator looks like 



2 ^ [COS(^) - COSh(^)]^ 2 , o2n 
(Jj 



(31) 



(32) 



and the position-dependent mass will be given by M {u, v) = mo a^/[cos(i') — 
cosh(M)]^. In this last example of anisotropic system, the effective potential 
happens to be 



V{u,v) 



a^mo 



a sinh(-u) 



cosh(K) — cos(v) 



a sm f 



cosh(ii) — cos(i') 



(33) 
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and the wavefunctions are correspondingly given by 

1 



^2(™+"-i)m!n!PiPs 



:H„ 



\/2 / a sinh('u) 
Pi \cosh(-u) — cos(t') 



X 



Hr, 



exp 



a sm(vi 



'V2 f 
P2 \cosh('u) — cos(f) 



exp 



a sinh(M) 



Pi \cosh(M) — cos(f ) 



a sm f 



P2 \cosh(M) — cos(f 
whose typical behavior appears in the Figure 4. 



4 The case of interaction with a magnetic 
field 

In this section, we first construct the SU{2) coherent state for the case of a 
particle with constant mass in the presence of a homogeneous magnetic field 
which, as far as we know, was not yet obtained in the literature. Then, we 
proceed by using that result in order to achieve the SU{2) coherent states for 
the PDM particles under the presence of some non-homogeneous magnetic 
fields. All the results will be obtained in the Coulomb gauge. Starting with 
the Hamiltonian 



H=-^t'--^(-t.f + t.-t) + ^t' + ^moCo' {x'+ y'), (35) 
2 niQ 2 tuq \ / 2 rriQ 2 

where e is the electrical charge and ~^{x,y) is the vector potential. 

Here we choose to work with a uniform magnetic field along the direction 
z, which is written as ^ = Pg ^- In the so-called symmetrical gauge it is 
written as 

'A = ^{-yi + xj), (36) 

Thus, beginning with the classical Hamiltonian defined in ( l35l) and quantizing 
it, one arrives at the following Schroedinger equation 



2 rriQ 



V2^ + 







1 h P, 

2 mo 



[X 



dy-y d,)^ + -mon^x^+ y^)^ = Etlj. (37) 
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with n = ujl + i^. 

U 8 mo 

In order to deal with the case of a constant mass in the presence of a 
homogeneous magnetic field, one can follow two alternative routes. One 
way is to write the system in polar coordinates and the another is keeping 
the cartesian coordinates and performing some convenient time-dependent 
transformations. This last is the one we will follow here [38] • For this, we 
can start with the time-dependent Schroedinger equation 



-VV + — ^ (x dy - ydx) a + U,ff a = th—, (38) 



2 mo 2 mo mg dt 

such that a = e~^^^'^'^^^xi^j u)- Then, we perform the time-dependent rotation 



Xi \ _ f cosa{t) sena{t) \ ( x 
Yi J y — sena{t) cosa{t) J \ y 

which, after some manipulations [M], and choosing a to be given by 



(39) 



where c is an arbitrary integration constant. One can finish with an effective 
isotropic two-dimensional harmonic oscillator such that 



2 mo 



+ Irno c^) {X^, + Y^) 



mo 2 



^ = '^%- (41) 



This allow us to map the original problem into one where the differential 
equation is given by 

-i^^'^'x + U.fix^Ex, (42) 

where x = x{Xi, Yi). Finally, using the usual variables separation procedure, 
one can arrive at the expression 



Xnn.{X,,Y^) = 0(Xi)^(ri) = 

= (j)[cos{a) X + sen{a)y] ip [—sen{a) x + cos{a)y] 
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(M [cos{a) X + sen{a) y, —sen{a) x + cos{a) y]) 



-1 



R V2"^+"-%n!m! 
y/2[cos{a) x + sen{a) y] \ 



X 



X H, 





(43) 



X H, 



X exp 



Finally, by using the above wave-functions, the SU{2) coherent states 
can be easily written. Their general aspect happens to be the same of the 
isotropic harmonic oscillator as presented in the left plot of Figure 2. 

Now, we can study the cases of PDM particles in the presence of magnetic 
fields. Following the same lines developed in the previous section, we make 
the variables transformation x = f{u,v) and y = g{u,v) in the Schroedinger 
equation and, since we already now the corresponding transformation in the 
Laplace operator, we only need to verify the impact of the transformation 
over the linear differential operator coming from which comes to be 

xdy-yd^ = j{[ f{u,v) fu + g{u,v) Qu] dy-[ f{u,v)fy + g{u,v)g^ ] du} . 

(44) 

Then, we finish with the transformed Schroedinger equation in the spatial 
variables u and v appearing as 



2M{u,v) 



2 M{u,v) 



iHBo 





(45) 



The magnetic field will be then given by 




(46) 
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Again, the first example to be analyzed is the one defined in Equation 
(IT^ . In this case we have that 

R{u,v) = [cl + cf) [v^ + v) - {Cidi + C2d2) u - {cid2 - C2di) V , 

S{u,v) = ^[cl + cl) [u^ + uv"^) + {c2di - Cid2)u + {cidi + C2d2)v. (47) 

Here, we will have a particle with mass M {u,v) = mo(cf + C2){u^ + f^) 
moving in the presence of an axially symmetric field with a quadratically 
growing intensity, 

t = Bo icl + cl){u^ + v^) t. (48) 

Choosing the parameters uiq = = ^2 = Ci = and C2 = 1, we restrict 
ourselves to the case of parabolic cylinder coordinates where 



M(u, v)=u^ + v'^ , R{u, v) = -^v {u^ + v^), S{u, v) = {u^ + v^), (49) 
and the potential governing the system is 



V{u,v) 



mo u) 



t[u^-v^Y + {u v) 



+ 



(50) 



and the wavefunctions are 

1 



R a/2™+"-i 



Hr, 



nnimi 



R 



R 



(51) 



xexp 



2 2\ 2 



2 R 



U V 

~R 



which, after plotting the probability density of the corresponding SU(2) co- 
herent state, presents the profile which is very similar to the one appearing 
in the Figure 1. 
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Now, in the case of elliptical cylindrical coordinates one have 

R{v) = cos(f) sin(t>), S{u) = cosh(M) sinh(u). (52) 
and the magnetic field is given by 

^ = a'^[cosh{2u) - cos(2t;)] 't . (53) 
The potential under which the charged particle is moving is 

Tfl 0? Up' 

V {u, v) = [sinh (iP) sin (t>)^ + cosh (m^) cos (f )^] + 

(54) 



■z T-, V [cosh(4 u) — cos(4 v) 
—e Bq 



2 [cos(2w) - cosh(2 u)] 

and the wavefunctions can be straightforwardly obtained through direct sub- 
stitution of the transformation functions in the expression of the wavefunc- 
tions f H3|) . After that, one can plot the probability density, which present 
the same general behavior as appearing in the Figure 3. There, one can see 
a sequence of regions where the highest probability density is over a closed 
curve, which refiects the periodicity of the transformation functions. 

In the third and last example, we will deal with the bipolar coordinates. 
In this case the Laplace operator is written as V^^ = (^^2_|_^2^^ 
and we get 

, cosh(u) sin(f) cos(i;) sinh(u) 

R{u,v) = - — rr^' -^^'^ = -T — r~\ TT^- 

[cos(f j — cosh(ujJ^ [cos^yj — cosh(MjJ^ 

the corresponding magnetic field looks like 

[cos(f j — cosh(njJ"^ 



and, in this last example, the potential is written as 

a^mo (sinh(u)2 + sin{vf) c^Bq , — — — 

V{u,v) = — —J - — — Vcosh(n)2 - cos{v)^. (57) 

^ (cosh(M) — cos(f )) 4 mo 

Once more, after constructing the wave-function for the SU (2) coherent 
state, one can see that we arrive at a Figure which is very similar to the case 
(a) of the Figure 4 below. 
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5 Final remarks 



In this work we explored a method for generating exactly solvable position- 
dependent mass particle in the present of some potentials. This approach 
was recently introduced in the case of one-dimensional systems [1]. Here we 
see that the extension for higher dimensions is not trivial and allows one to 
analyze complex and interesting new systems. Moreover, we conduce our 
study by computing the so called SU (2) coherent states, which are specially 
interesting wave-packets which present their maximum of probability over 
the classical trajectory [H [391 SOI SI]- As particular examples of variable 
transformation, we used a polynomial example as well as the cases of the 
cylindrical elliptical and the bipolar coordinates. 

In the case of a charged particle under the effect of pure homogeneous 
external magnetic field, despite the fact that the we are dealing with time- 
dependent wave-function solutions, it was verified that the corresponding 
SU{2) coherent states do not present dependence in the time variable. In 
that case, the maximum of the probability density is over a circle. Then, per- 
forming the change in the spatial variables, we got those states for the case of 
some particular spatially dependent masses and non-homogenous magnetic 
fields. In the case where the transformation used was that of polynomial 
form, it was observed that when the zero degree coefficient of the polynomial 
used in the transformation was non-zero, the maximum of probability splits 
into two ones which become more and more distant from each other when 
that parameter increases. Furthermore, in the case of the other two transfor- 
mations used here, the probabihty density was naturally periodic as a direct 
consequence of the periodicity of the transformation. 
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Figure 1: Probability density for the case of polynomial solutions. Here the 
parameters used are a — Ci — C2 — 1 and: a) d2 — ^ and c) d2 — 7. 




-4-20 2 4-4-20 2 4 



Figure 2: Probability density for the case of polynomial solutions. Here the 
parameters used are: di = d2 = c-i = 0, C2 = 1 and the frequencies are the 
following: a) cui — 1 and 002 — 1 and c) cui — 2 and 002 — S. 
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Figure 3: Probability density for the case of elliptical cylindrical coordinates. 
In this case the frequencies were: a) oui — 1 and 002 — 1- 
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Figure 4: Probability density for the case of bipolar coordinates. The used 
frequencies were: a) cui — 1 and 002 — 1, b) cui — 1 and 002 — 2 and c) oui — 2 
and UJ2 — 3. 
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